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Abstract 



We develop a shell-model Monte Carlo (SMMC) method to calculate den- 
^" , sities of states with varying exciton (particle- hole) number. We then apply 

\^ \ this method to the doubly closed-shell nucleus ^''Ca in a full Qs-ld-^f-lp 

^D ' shell-model space and compare our results to those found using approximate 

^— V ' analytic expressions for the partial densities. We find that the effective one- 

ON ■ body level density is reduced by approximately 22% when a residual two-body 

interaction is included in the shell model calculation. 
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1. Introduction 

Particle-hole, or exciton, level densities enter into the description of partial decay rates in 
nuclear preequilibrium emission ll|-0. These level densities have been modeled using analytic 
expressions [^,Q that describe nuclear excitations in terms of the number of particles, p, and 
holes, h, measured from the Fermi surface, with the exciton number Ne = {p + h)/2. 

For a single species of particles, Williams |^ derived an expression for the partial density 
of states given by 

where E is the excitation energy measured above the ground-state configuration, g is the 
single-particle density of states, and G/g plays the role of an effective Pauli energy with 
G = (p^ + /i^)/4+(p — /i)/4 + /i/2. There exist more complicated expressions that distinguish 
between protons (with single-particle density gp) and neutrons (gn) in a given nucleus, but 
we will not quote them here. In the most naive picture, a uniform spacing of single-particle 
states, d, is assumed, in which case the single-particle level density is g = 1/d, measured in 
units of MeV-^ 

Equation |I| and its neutron/proton counterpart suffer from several deficiencies including 
the assumption of an unlimited number of single-particle states, an inexact treatment of the 
Pauli principle, and the assumption of a uniform single-particle level spacing. Extensions to 
the basic model that ameliorate some of these effects have also been pursued. As examples, 
we mention Bogilia et al. in which the energy dependence of the single-particle level 
spacing was included in a general way; using the equidistant single-particle picture, Kalbach 
and Zhang and Yang considered Pauli principle corrections to the state densities; and 
De and Hua [Q considered the effects of pairing in addition to the Pauli blocking on the 
state densities. These effects were combined and extended to non-uniform level spacings by, 
for example, Harangozo et al. [l^ . 



A further difficulty with the simple formula is that the residual two-body interaction, 
which is present beyond the nuclear mean field, and which includes important contributions 
beyond J = pairing, is not incorporated. While some progress has been made to approxi- 
mate the effects of the residual interaction |Tl| on the partial level densities, no interacting 
shell-model calculations have been performed in large model spaces that would indicate the 
effect of the two-body interaction on the single-body density parameter, g, nor have there 
been any partial density-of-state calculations in the interacting shell model. 

In this article, we describe calculations that study the effects of the residual two-body 
interaction on the partial level densities, and present results for partial level densities in 
"^^Ca. Our approach is to study a related quantity, Yn^{P), which is the ratio of the particle- 
hole partition functions, Zn^{P), to the full partition function, Za{P), as a function of the 
inverse temperature, f3, (measured in MeV~^) in the system. We perform our calculations in 
a full Os-ld-Of-lp model space using shell-model Monte Carlo (SMMC) techniques [|12|,|13 



and an interaction that describes reasonably well the low-lying spectral properties of nuclei 
in the sd-fp region [|14|. We will compare our results with those obtained from eq. (|lD and 



its proton/neutron counterparts. Finally, we will show partial densities of states for several 
exciton numbers. In section 2, we give an overview of our calculational method. We present 
results in section 3, and conclude with a brief summary in section 4. 



2. Calculation of Excitons in SMMC 

Investigations into both ground-state and thermal properties of nuclei have been de- 
scribed using the SMMC technique |]12|. This method offers an alternative description of 



nuclear structure properties in the shell-model context that is complementary to direct diago- 
nalization. SMMC is designed to give thermal or ground-state expectation values for various 
one- and two-body observables. Indeed, for larger nuclei, SMMC may be the only way to 
obtain information on the thermal properties of the system from a shell-model perspective. 
In this method, we make use of the imaginary time many-body propagator tl = exp{—(3H) 
to calculate the expectation values. For example, the excitation energy of a nucleus is 
^iP) = {H{P)) — {H{oo)), where {H{oo)) is the ground-state energy. In order to find the 
excitation energy of a nucleus with particle number A, we must then calculate 

TtPaU TtaU 

where Pa = S{N — A) projects the trace over all many-body states in the system to those 
states that have the desired particle number. 

Two-body terms in H are linearized through the Hubbard- Stratonovich transformation, 
which introduces auxiliary fields over which one must integrate to obtain physical answers. 
Since H contains many two-body terms that do not commute, one must discretize (3 = NtAjS. 
The method can be summarized as 

r - iNt 

Za = TtaU = TTAexp{-(3H) -^ Tta [exp{-A(3H)\ 

Nt 



fv[a]Gia)TTAl[exp[A(3h{a^)] , (3) 

•' n=l 



where o"„ are the auxiliary fields at a given imaginary time-step A(3 (there is one cr-field 
for each two-body matrix-element in H when the two-body terms are recast in quadratic 
form), V[a] is the measure of the integrand, G{a) is a Gaussian in a, and h is a. one-body 
Hamiltonian. Thus, the shell-model problem is transformed from the diagonalization of a 
large matrix to one of large dimensional quadrature. Dimensions of the integral can reach up 
to 10^ for systems of interest, and it is thus natural to use Metropolis random walk methods 
to sample the space. Such integration is most efficiently performed on massively parallel 
computers. Further details are discussed in Koonin et al. ||12|| . 

In order to obtain density-of-state information, we calculate in SMMC the expectation 
of the energy and integrate the thermodynamic relationship 

Em-'-^ (4) 



to obtain 



r/3 



In Za{(3) = - f d(3'E{(3') - In Za{0) , (5) 

Jo 

where Za{0) = Tr^l is the total number of A-particle states in the system. Za{P) and p{E) 
are related by the inverse Laplace transform 



/oo 
dEeM-l3E)p{E) , (6) 

-oo 



which can be solved in a saddle point approximation to yield 

exp(S') 



p{E) 



V27r/?-^C 



S = ^E + \nZ{^)- I3-'C = -—. (7) 

In this expression C is the heat capacity of the system. 

For our discussion, we study ^°Ca. The non-interacting ground state is a filled sd- 
shell with no particles in the fp-shell. Our excitons are then enumerated with respect 
to the filled sd-shell. We may excite both protons (tt) and neutrons (u) so that A^e = 
(pvr + hn + Pu + hy)/2 = Nfp, where Nfp (Nsd) gives the number of particles in the fp- (sd-) 
shell. For example, Ne = 2 includes the following particle- hole excitations: {OpnOh-,,, 2py2hv) 
(Ip^l/i^, Ipylhy) {2pT,2h^, Op^Ohi,). Furthermore, < A^g < 24 since, at most, 24 particles 
can be excited from the sd-shell into the /p-shell. 

The ratio of the partition function for N^. excitons to the full partition function for the 
^-particle system may be found by introducing a second number projection operator, 

Pn. = S{Nsd - Nsd)S{NfP - Nfp) , (8) 

provided that A = Nsd + Nfp. In reality we perform this projection for both protons and 
neutrons simultaneously, but this only complicates notation and will not be discussed here. 
We calculate the ratio of partition functions as 

Therefore, J^n^ ^Afe(/^) = 1 which we use as a convenient numerical check. We also extract 
the energy of the particle-hole excitations, E]\f^ , as 

E,,W^y-}^ + EiP). (10) 

We may now employ eq. (0) for the partial density of states: 

exp(S'ArJ 



PN^E) 
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S^^=P^^E^^+lnZ^M; (^NlCNe = -'^. (11) 

Here /^tv^ = PNeiE^J is determined by inverting the relation E^^ = EN^iPNe): and Cn^ is 
the heat capacity for the particular exciton number. 



3. Results 

We now turn to a description of our ^°Ca calculation in the Os-ld-Of-lp shell model 
space. Our starting point for an appropriate interaction is taken from ref. [|l^]. In order to 
obtain a microscopic effective interaction, we begin with a free nucleon-nucleon interaction 
which is appropriate for a description of low-energy nuclear structure. The choice made 
in ref. |T^ was to work with the charge-dependent version of the Bonn nucleon-nucleon 



potential model as found in ref. |]I5[. Standard perturbation techniques, as discussed in 
ref. |]TB[, were employed to obtain an effective interaction in the full sd-fp model space. 
Finally, the interaction was modified in the monopole terms using techniques developed by 
Zuker and co-workers JT7|p!B| . 

SMMC calculations for realistic interactions typically have a Monte Carlo sign problem 
which can be overcome by an extrapolation technique discussed in ref. |l^, and successfully 



applied to the sd-fp region in [Q. This extrapolation technique was also applied to thermal 
properties of nuclei |^0|], but the statistical error inherent in the energy upon extrapolation 
prevents a full description of the density of states unless one has good justification to spend 
the computational resources to reduce the statistical error. It was recently demonstrated 
that a good reproduction of the experimental density of states could be obtained for nuclei 
in the Oflp-OgQ/2 shell |21[] using a pairing-plus-quadrupole interaction that was free from 
the sign problem. In this work, we fit our realistic two-body interaction discussed above to 
a pairing-plus-multipole interaction given by 

H2 = gonP^,Poo + ^nJ2xu: Y.^-YQ'^.Qu-, : , (12) 

where :: denotes normal ordering and P|^, Q,yf^ are pair and quadrupole operators given by 

pi=j:i-y''Ua\\y.\\jb)[dixdi],, 

ab 

In eq. (|13D a = nij denotes a single-particle orbit and djm = {~y~^"^dj-m- We fit go and 
the /i = 2,4,6 multipoles to the realistic interaction from ^^. A least-squares fit gives 
an interaction which indeed has a good Monte Carlo sign. After some minor adjustments 
to the pairing strength in order to obtain a better gap between ground states and first 
excited states in several light nuclei, we use the following parameter set: go = —0.63 MeV, 
X2 = -0.047 MeV fm-^ X4 = -0.001 MeV fm'^ and Xe = -0.17 x lO'^ MeV fm'^^^ 
(Large enhancements of the {a \ r'^ \ h) matrix elements as v increases is the reason for 
the decrease in Xv values, although contributions to two-body matrix elements arising from 
the higher multipoles is significant.) Our single-particle energies are 0.0, 5.36, 0.64, 8.21, 
14.21, 10.14, and 12.07 MeV for the 0^5/2, 04/2, lsi/2, O/7/2, O/5/2, lp3/2, and lpi/2 orbitals, 
respectively. We do not correct for center-of-mass motion in these calculations, although 
such a contamination to the Y^^ should be fairly small in this system. Furthermore, we 
do not include odd multipoles which upon fitting were found to give coefficients that cause 
Monte Carlo sign problems. Thus our negative parity states are probably less well described 
by this choice of Hamiltonian. 



We use the results of the non-interacting case to demonstrate the vahdity of our technique 
for finding the partial partition functions. In order to show this, we calculate by enumeration 
the total number of many-body states for each N^. This can most easily be done by using 
eq. (ID to find Za{P = 0). We also find Yjy^^O) by an extrapolation from small, but finite, j3. 
We show in fig. |1] our results for the number of states as a function of the exciton number. 
The SMMC results are compared to an exact counting of the number of states of a given 
exciton number. The agreement is excellent. The total number of calculated SMMC states 
is 3.834x10^^ as compared to the exact value of 5.095x10^^. 

We show in fig. ^ a comparison of the non-interacting (left) and interacting (right) 
calculation. The A'^e = calculation gives some indication of the thermal freeze-out of 
the ground state. The non-interacting calculation requires fairly large /3 to fully reach 
the ground state, since the first excited state is only 0.64 MeV above the ground state. 
We pursued these calculations to /3 = 4.0 MeV~^, for which (H) is 0.139 MeV from the 
ground-state value. Since it takes more thermal energy to overcome the pairing interaction 
and to excite them from the ground-state configuration, the interacting curves representing 
different Ne excitations are compressed in jS relative to the non-interacting curves. (The 
excitation energy of the first excited state is approximately 3.5 MeV.) However, the same 
features remain. Clearly, as thermal energy is decreased, and f5 becomes larger (temperature 
decreases), it is more difficult to produce large particle- hole excitations in the system. The 
converse is also true at higher temperatures, where it is difficult to obtain only N^, = 2, for 
example. In the interacting case, the low temperature tail of the Ne = 2 exciton tends to 
spread further in (3 than does the Ne = 1 tail. As we shall see, this has direct consequences 
on the partial densities of states. Furthermore, since J^n^ '^NeiP) = 1; "we can interpret Y^^ 
as a measure of likelihood to find Ne excitons at a given temperature. Since the excitation 
energy is a monotonic function of the temperature, one expects the density of states to be 
dominated by excitons of a particular type in a given energy range. As we shall see, this is 
indeed the case. 

We compare our results to those obtained from eq. (0) and its proton/neutron equivalent 
by adjusting g to obtain a fit to the calculated Y^^ curves. This is demonstrated for the 
A'e = 4 curve in fig. ^. We fit each of our SMMC curves to eq. (|I|) for an effective one- 
body density parameter g^s = gp + gn [HI, and the effective level density parameter is 



a = (7r^/6)5'cff- A uniform Fermi gas yields a ~ A/15 (=2.67 for A=40) MeV~^, a harmonic 
oscillator potential yields a ~ A/10 MeV~^ (=4.0), and the empirical value is A/8 MeV~^. 
We obtain a = 5.43 {ges = 3.3) MeV~^ for the non-interacting case, and a = 4.27 {ges = 2.6) 
MeV^^ in the interacting case, rather independent (within 0.02 MeV~^) of A^e > 2. For 
A^e = 1 the comparison cannot be made as the Blann- Williams formula breaks down. Thus, 
ges is reduced by 22% in the presence of an interaction. 

The decomposition of the A"e = 4 case into the various proton-neutron components is also 
shown in fig. ^. The (Op^O/i^r, Ap^Ah^) plus its neutron counterpart carries very little weight 
here, while the (2p^2/i^, 2p^2h^) component of Y^ carries the most weight. As expected, in 
all Ne cases the largest component of Y/v^ is the one in which the number of excited neutrons 
equals the number of excited protons. 

Finally we present our result for the calculation of p^^ derived by using Z^^ (/3) in eq. (^ 
and En^ from eq. (p^O]). The natural log of pn^ is shown in fig. H for the A'e=1^6 excitations. 
We also include in the figure the total state density as a function of the excitation energy 



in the system. The saddle point approximation breaks down in regions where there are 
very few states, which makes it difficult to describe well-separated states in the low-lying 
spectrum {E* < 3 MeV) for the full density or for the individual exciton densities. We also 
propagated our statistical error bars through the calculation of p^^, but, as can be seen, 
they are very small except in the case of the A'"e = 1 excitons. Note that the majority of 
states, for example at E* = 25 MeV, are A^e = 4 states, while at E* = 35 MeV the A^e = 5 
states contribute most. This localization in excitation energy of excitons was reflected in our 
earlier discussion of the behavior of Yat^. Interestingly, the N^, = 1 density of states begins in 
energy slightly above 2p2h state density. Recall that experimentally the ffist excited state of 
"^°Ca is a 0^ 2p2h state (at 3.3 MeV), and that the ffist negative parity state (a 3^) occurs at 
a slightly higher energy of 3.7 MeV. Our Hamiltonian fairly closely gives the correct relative 
starting energies for these two exciton configurations, although, due to the breakdown of 
the saddle-point approximation for low state densities, we cannot precisely determine the 
excitation energy of the ffist excited O"*" level. We also note an interesting pairing effect 
that shows up in the partial densities. Note that the N^ = I state density starts about 
0.6 MeV above the N^ = 2 case. The A'^e = 3 state density begins approximately 2 MeV 
above the iVe = 4 case. This is a manifestation of pairing in the system. It takes more 
excitation energy to produce an odd particle-hole excitation than it does to produce an even 
particle-hole excitation since energy must be expended to break pairs for odd excitations. 
This effect was already apparent from the discussion of the low-temperature behavior of the 
Yat^, as indicated above. 

4. Conclusion 

We have discussed in this article how one may obtain information on particle-hole ex- 
citations using SMMC methods for calculations of a nuclear system. Our technique uses 
the ratio of the particle-hole partition function to the full partition function of the system. 
The method incorporates exact Pauli blocking, non-equidistant single-particle energies, and 
gives the exact partial densities for a given nuclear effective interaction, within statistical 
errors. It also has a well-defined energy scale. One drawback of the present calculation is 
that the space size is limited to two major oscillator shells, although this can be rather easily 
overcome. The projection operator introduced in this may be applied in any Monte Carlo 
technique where ratios of partition functions are needed. Our results also indicate that the 
effective (/-parameter used in eq. (1) should be reduced by approximately 22% to account 
for the inclusion of the two-body interaction which acts to correlate the nucleus beyond the 
simple mean-field or pairing prescription. 

The method we have described here could be further advanced in two ways. One may 
increase the model space used, thus allowing for a broader range of energies and excitation 
modes to be explored. The method is also applicable to mp-nh excitations if we extend our 
studies to open-shell nuclei such as, e.g., ^^Ca. This may be pursued in future work. 
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FIGURES 

FIG. 1. Top panel: Running sum of the non- interacting states found from Z^r^ at /3 = 
compared with the running sum from an exact calculation. Bottom panel: the total number of 
states for each N^. found in both non-interacting SMMC calculation and compared to the exact 
result. 

FIG. 2. A comparison of the interacting (right) and non-interacting (left) functions Yat^ for 
A^e = 0, • • • , 6. 

FIG. 3. Decomposition of the Yjs[^=4, into its particle hole components, and a comparison of 
the full yivg=4 generated from SMMC to that obtained from eq. (||). The fit corresponds to 
5eff = 2.6MeV-i. 

FIG. 4. Calculated partial densities of states p^S^Ne) for N^. = l,---,6, using the saddle 
point approximation from eq. (0). Also plotted is the total state density calculated from eq. ^. 
Statistical error bars that are not visible are smaller than the symbols used. 
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